Site-Dilution in quasi one-dimensional antiferromagnet Sr2(Cui_ x Pd x )Os: reduction of 
Neel Temperature and spatial distribution of ordered moment sizes 
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We investigate the Neel temperature of Sr2Cu03 as a function of the site dilution at the Cu 
(S — 1/2) sites with Pd (S = 0), utilizing the muon spin relaxation (/iSR) technique. The Neel 
temperature, which is Tjv = 5.4K for the undoped system, becomes significantly reduced for less than 
one percent of doping Pd, giving a support for the previous proposal for the good one-dimensionality. 
The Pd concentration dependence of the Neel temperature is compared with a recent theoretical 
study (S. Eggert, I. Affleck and M.D.P. Horton, Phys. Rev. Lett. 89, 47202 (2002)) of weakly 
coupled one-dimensional antiferromagnetic chains of S = 1/2 spins, and a quantitative agreement 
is found. The inhomogeneity of the ordered moment sizes is characterized by the /J.SK time spectra. 
We propose a model that the ordered moment size recovers away from the dopant S = sites with 
a recovery length of £ « 150 — 200 sites. The origin of the finite recovery length £ for the gapless 
S — 1/2 antiferromagnetic chain is compared to the estimate based on the effective staggered 
magnetic field from the neighboring chains. 

PACS numbers: PACS numbers: 76.75.+i, 75.25.+Z, 75.10.Jm 



I. INTRODUCTION 

The discovery of high-T c cuprates has prompted 
theoretical and experimental investigations of low- 
dimensional spin systems with spin quantum number 
S = 1/2. There were a series of neutron diffraction 
studies reported in the high-T c cuprate La2-a;Sr x Cu04 
0, , in which the dynamical spin correlation was mea- 
sured. It was found that the inverse correlation length 
^ 1 (x,T) may be divided into the two terms: the tem- 
perature independent term £ -1 (a;,0) which is purely de- 
termined by the doping concentration and the doping 
independent term £ _1 (0,T) which follows the universal 
temperature dependence. In samples with Neel ordering, 
the correlation length diverges £ _1 — > at T/v, exhibiting 
long-range magnetic order. Recently, effect of static site 
dilution was also investigated in La2(Cui_ x (Mg,Zn) x )04 
[j| . The Neel temperature was found to disappear at the 
classical percolation threshold p c — 0.407, but the site- 
dilution dependence of T/v does not follow the mean-field 
calculation. The reduction of T/v, spin-stiffness p s and 
the equal-time correlation length £(x, T) were compared 
with the microscopic quantum mechanical calculations of 
the S =1/2 Heisenberg model on a square lattice Q. 

The localized holes in the lightly doped 
La2- K Sr 2; Cu04 and the impurities in the site dilu- 
tion in La2(Cui_2,(Mg,Zn)2,)04 break the translational 
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symmetry. The ordered moment size of these systems 
should not be homogeneous in space. However, the 
signature of the spatial inhomogeneity of the ordered 
moment sizes in cuprates is difficult to detect in the 
momentum space, because the doped sites are randomly 
distributed, the spin quantum number S is small and 
the dimensionality is low. Consequently, the length scale 
relevant to the moment size distribution in the impurity 
doped cuprates has not been resolved by the previous 
neutron diffraction measurements. 

Site dilution in the quasi one-dimensional spin sys- 
tems exhibits more important features of the quantum 
spin systems. One peculiar effect of site dilutions in one- 
dimensional spin systems is the creation of a Neel order 
out of the singlet ground state, as discovered in the im- 
purity doped spin-Peierls material CuGeO-3 0,0 and the 
2- leg spin ladder material SrCu203 @. These discover- 
ies have promoted the idea that the Neel state appears 
as a competing phase to the original spin-gapped singlet 
state of either spin-Peierls or the spin-ladder. As a result 
of competition, there exists a length scale which deter- 
mines the spatial variation of the ordered moment sizes 
||. In doped spin Peierls compounds, the spin de- 
cay length in the Neel ordered state is determined by the 
ratio between the spin-gap magnitude A and the intra- 
chain antiferromagnetic coupling J, such as £/<x = J/ A 
• This suggests that the spatial inhomogeneity of the 
ordered moment size is a feature of the spin-gap, and may 
not appear in the regular quasi one-dimensional antifer- 
romagnets with Neel order, which are gapless (A = 0) 
and the correlation length diverges. 
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Recently, a series of measurements of Cu benzoate, a 
compound with 5 = 1/2 one-dimensional chains, have 
identified the existence of an energy gap induced by the 
external magnetic field 

UGH 

The origin of the energy 
gap has been interpreted as the effect of the staggered 
magnetic field induced on the Cu sites by the combi- 
nation of the anisotropic g-tensor and the external field 
[llj . This observation points out that the 5 = 1/2 spin 
chain, which has a singlet ground state and gapless exci- 
tations if isolated, may acquire a spin-gap under a certain 
type of perturbation. In quasi one-dimensional antiferro- 
magnets, the most common perturbation to an isolated 
chain is the existence of inter-chain interactions J'. In the 
mean-field approximation, the inter-chain interaction in- 
duces an effective staggered magnetic field B as zJ'(s z ), 
where (s z ) and z are the magnetic order parameter and 
the number of the nearest neighboring chains, respec- 
tively. The inter-chain interaction causes the Neel order 
[l2^ which is gapless because of the translational sym- 
metry of the system. With the impurity doping to the 
5=1/2 spin-chain, the translational symmetry is bro- 
ken and the hidden features of the effective staggered 
fields from the inter-chain interaction may appear in the 
form of the recovery length £ of the ordered moment size. 
However, this problem is still an open question. 

The site dilution should have a destructive effect in 
the long-range coherence of the Neel order, especially in 
low dimensions. The true one-dimensional chain becomes 
fragments of finite sizes upon impurity doping, and the 
long-range order becomes impossible. However, with the 
existence of inter-chain interactions, the destruction of 
magnetic ordering with site dilution may be moderated; 
the true disappearance of the Neel order may occur only 
at the percolation threshold for the three-dimensional lat- 
tice structure of the inter-chain interactions. Recently, 
Eggert et at presented a theoretical estimate for the 
Neel temperature as a function of the site dilution in 
the 5=1/2 quasi one-dimensional antiferromagnet. The 
predicted Tjv exhibits a simple reduction as a function 
of the average chain length L. Eventhough the ground 
state of an isolated 5=1/2 chain is a singlet, there was 
no enhancement of the Neel temperature upon impurity 
doping, in contrast with the gapped 5=1/2 systems, 
such as the spin-Peierls systems and the spin-ladders. 

In order to experimentally investigate the theoreti- 
cal predictions about Jjv, and the uniformity of the or- 
dered moment sizes in the depleted spin-chains, we have 
performed an investigation of Pd-doped Sr 2 Cu0 3 . The 
cuprate Sr2Cu03 has received an attention as a model 
material of the 5=1/2 quasi one-dimensional antifer- 
romagnet; its low Neel temperature Tjv = 5.4K 03 and 
the large in-chain interaction J ~ 2200K [D^ llfl Irn j sug- 
gests its good one-dimensionality. The ordered moment 
size « 0.06^s has been obtained by neutron scattering 
and /iSR measurements |14| . This value is strongly re- 
duced from the full moment size (= I^b) expected for the 
5 = 1/2 spins. The suppressed moment size of Sr 2 Cu03 
follows the prediction based on the "chain mean-field" 



theory 12], which employs the rigorous results of the iso- 
lated chains and includes the weak inter-chain interaction 
as the mean-field. As the non-magnetic impurity at the 
Cu site, Zn or Mg substitutions are the first choice, as has 
been already performed in the two-dimensional cuprates 
However in S^CuOs, Zn or Mg ions does not go into 
the Cu site. We have employed Pd ion instead; S^CuOa 
has an isostructural compound Sr2Pd03, in which the 
Pd 2+ ions are in the low-spin (5 = 0) state [H| • This 
compound has enables us to investigate the non-magnetic 
impurity doping to the 5=1/2 antiferromagnetic spin- 
chain. 

The structure of this paper is as follows. In section 
UTI we present the magnetic susceptibility data and the 
result of muon spin relaxation (/zSR) measurement of 
Sr2(Cui_2;Pd x )03. The Neel temperature was estimated 
from the temperature dependence of the muon relaxation 
rate. In section ITTT1 we calculate the magnetic field dis- 
tribution expected for the depleted spin-chains with Neel 
order. We assume a zero-moment at the impurity site 
and a recovery length £ to describe the recovery of the 
ordered moment size into the bulk chain. The calcu- 
lated field distribution is Fourier transformed to obtain 
the /xSR spectrum, and employed in the analysis of the 
nominally pure as well as the doped samples. Conclu- 
sions are presented in Section llVI 



II. EXPERIMENTAL RESULTS 

We grew single crystals of Sr 2 (Cui_ :E Pd 2 ;)03, employ- 
ing the traveling-solvent floating-zone technique, with 
CuO as the solvent. Stoichiometric ratio of SrC03, CuO 
and PdO powders are prepared, mixed in a mortar for an 
hour and pre-fired in air at 900°C in an AI2O3 crucible. 
The powder sample is again mixed and fired in air at 
950°C, before formed into a pressed rod using a rubber 
tube and a water static press. The polycrystalline rod is 
fired on a Pt plate at 1050°C. It is important to harden 
the rod by firing at the highest temperature possible for 
a stable growth of single crystal in the floating-zone fur- 
nace. We employed a gold-mirror bi-focus furnace made 
by NEC. 

The magnetic susceptibility of Si^Cui-zPd^Os crys- 
tal is shown in Fig. with the magnetic field applied 
parallel to the longest crystallographic axis. An increase 
of the Curie- Weiss component was observed as the Pd 
concentration increases. We assume the conventional 
form for the magnetic susceptibility: 



X(T) 



C 



r + e 



ir 



Xo 



(1) 



where C is the Curie term, ®\y is the Weiss temperature 
and Xo describes the temperature independent suscepti- 
bility as a sum of Van Vleck paramagnetism and core 
diamagnetism. 

The parameters to describe susceptibility are shown in 
Table HJ The Curie term is shown as the concentration of 
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FIG. 1: (a) Magnetic susceptibility of single crystalline 
Sr2(Cui_ I Pd I )03, with the magnetic field applied parallel 
to the longest crystallographic axis, (b) Pd concentration de- 
pendence of the Curie term. We note that the substitution 
with one Pd ion creates half an impurity moment assumming 
tht the induced moment size is \\xb- 



TABLE I: Parameters of susceptibility and /^SR 
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"obtained from an analysis of /iSR spectrum (see Discussion), 
''numbers in parenthesis are estimates from the Curie term. 



the impurity moments which are assumed to be lfis- The 
calculated concentrations of impurity moments are about 
half of the doped Pd concentration, as shown in Fig.^i. 
This Pd concentration dependence of the Curie term may 
be explained as follows: One Pd ion creates one chain 
fragment. Assuming that the intra-chain interaction J 
is much larger than the interaction between the chain 
fragments, the total spin quantum number of the chain 



fragment is well defined: it takes either S — or S = 1/2 
value, depending on the length of the fragment being 
an even or odd number of spin sites, respectively 
The observed Pd concentration dependence of the Curie 
term is consistent with the idea that it originates from 
the total spin of the created chain fragments. The Weiss 
temperature Qw exhibits a non-monotonic dependence 
on the Pd concentration. The origin of this dependence 
is unknown, however, its temperature scale is at most 
Qw ~2K, which is being significantly smaller than the 
intra-chain interaction J ~ 2200K. This suggests that the 
interaction between the chain fragments are negligible 
compared to the intra-chain interaction, satisfying the 
condition assumed in the theory ^| , which is compared 
with our measurement in the discussion section. 

It is known that the Neel order of S^CuOa is not de- 
tectable by magnetic susceptibility ^3] • This is probably 
because of the very small ordered moment size w 0.06/ib. 
This feature of the material requires employing the muon 
spin relaxation (/iSR) technique, which has the highest 
sensitivity among other experimental techniques to de- 
tect the magnetic order with small and/or dilute mag- 
netic moments. We performed zero-field muon spin relax- 
ation measurement on Sr2(Cui_ 2; Pd a ;)03 crystals at M15 
beam line of TRIUMF (Vancouver, Canada). Muons 
with 100% spin polarization were injected into the single 
crystalline sample with the initial polarization parallel to 
the longest crystallographic axis. This geometry is the 
same as the one employed in the previous measurement 
of nominally pure Sr2Cu03 14]. The time evolution of 
muon spin polarization in zero-field is shown in Fig. 

In undoped Sr2Cu03, muon spin precession was ob- 
served as a consequence of Neel ordering 14] . This in- 
dicates that the ordered moment size is relatively homo- 
geneous in space, and the local field at the muon sites 
is well defined. Upon Pd doping, the muon spin preces- 
sion disappears and the time evolution of muon spin is 
dominated by the relaxation which has an approximately 
exponential behavior as a function of the time. Such ex- 
ponential relaxation signal in zero-field comes from the 
1/Ti relaxation caused by spin fluctuations, or alterna- 
tively from static fields of spatially distributed magnetic 
moments. The dynamic and static situations for the 
muon spin relaxation may be distinguished by the "de- 
coupling" measurements under longitudinal fields applied 
parallel to the initial muon spin polarization[19j. The re- 
sults are shown in Fig. |2Jd for the x = 1.0% specimen, 
together with the analysis us ing the static relaxation in 
Lorentzian field distribution 20] . It is clear that the re- 
laxation is caused by the static field distribution, as is 
evident from the time-independent behavior in the long 
terms [20]. This leads to the conclusion that the muon 
spin relaxation is caused by a magnetic order. 

We analyze the zero- field muon spin relaxation P M (t) in 
the doped samples utilizing the phenomenological func- 
tion: 

P„(t) = exp(-i(A„ d i) 2 )exp(-A mQff i) (2) 
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FIG. 2: (a) Zero-field muon spin relaxation in 

Sr2(Cui_ I Pd I )03. The solid lines for the doped samples are 
the fits using eq. @. On the spectrum of the samples with 
x = 0, the solid line is a plot of eq.(l) in Ref . . The dashed 
line is the sum of two exponential functions to describe the 
background relaxation, (b) Muon spin relaxation in longitu- 
dinal field (H || P M (t)) for the x = 1.0% sample. The solid 
lines are the Lorentzian functions in the longitudinal fields 
used in the analysis. 
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FIG. 3: Temperature dependence of the exponential relax- 
ation rate for S^Cui-zPdxOa. The Neel temperatures are 
estimated from the temperature points where the relaxation 
rate start to rise. 



as shown by the solid lines in Fig.^. In the analysis, the 
Gaussian relaxation was assumed to describe the static 
nuclear dipolar- fields, and was set to A n d = 0.12/xs -1 
independent to the temperature and the sample. The 
exponential relaxation rate A TOQff , which parametrizes 
the muon spin relaxation caused by the (atomic) mag- 
netism, are plotted in Fig.|3as a function of temperature. 
The Neel temperature was defined as the temperature at 
which the exponential relaxation rate start to increase; 
the estimated Neel temperatures are shown in the figure. 



III. DISCUSSION 

The absence of muon spin precession in Pd doped sam- 
ples indicates the fragility of the spatially homogeneous 
ordered moment sizes in Sr2Cu03. With only less than 
one percent doping of Pd, the coherent precession of 
muon spin disappears, indicating that the ordered mo- 
ment sizes have a broad spatial distribution under site 
dilution. The distribution of the moment sizes may be 
broader than in Zn- or Mg-doped CuGeOa, where muon 
spin precession was observed in the Neel state together 
with the exponential relaxation |2l|. In the doped spin- 
Peierls compound, the distribution of the local field was 
consistent with the model that maximum moment size 
is induced near the doped center, and the moment size 
reduces exponentially with the coherence length ks 10 
spin sites |2l|. In Sr2(Cu,Pd)03 where the Cu moments 
are depleted by Pd impurities, the maximum moment 
should be located in the middle of the chain fragments 
as shown in Fig. 0] Because of the long coherence length 
expected for the gapless 5=1/2 spin chain, the effect of 
the non-magnetic impurity may be extended in a large 
area; the maximum moment size of one chain-fragment 
might strongly depend on the chain length, and may not 
have a well-defined value in the doped system. 

In E.ef . [l4j . the /xSR signal of the nominally undoped 
Sr2Cu03 was analyzed by assuming two muon sites each 
of which consists of one precession signal and one expo- 
nentially decaying signal. The latter was interpreted as 
1/Ti relaxation of the local field component parallel to 
the initial muon spin orientation. The exponentially de- 
caying terms describe the background relaxation which 
exhibits a fast front-end before f» 1/Lts as shown by the 
dashed line in Fig. [5^. The corresponding relaxation rate 
for the front-end (« 5/xs -1 ), however, is too large for the 
residual dynamics in the Neel ordered state. In this sec- 
tion, we calculate muon spin relaxation for the inhomo- 
geneous Neel order with the existence of spin vacancies, 
from which the ordered moment size recovers exponen- 
tially with the recovery length £ as shown in Fig.0J The 
vacancies are the doped Pd ions or the impurities remain- 
ing in the nominally pure sample which manifest them- 
selves as the Curie term in the magnetic susceptibility 
(Fig. QJ. For the nominally pure Sr 2 Cu03, the Curie 
term corresponds to the impurity level « 0.04% which 
corresponds to L w 1/(2 x 0.04 x 10~ 2 ) « 1000 sites of 
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FIG. 4: Spatial evolution of the magnetic moment size as- 
sumed on a chain fragments. The small circles are the doped 
Pd sites, which were assumed to behave as an S — impurity. 



unperturbed spin-chain. 

With the model shown in Fig.0| the moment size S(z) 
behaves as a function of the distance z from the end for 
a chain fragment with length Lq: 



S(z) = ,s 1 - 



cosh((z - |£oVO \ 
3osh(±L /0 J 



(3) 



where sq is the moment size at infinity (sq ~ 0.06/is) 
and £ is the moment-size recovery length. By assum- 
ing that the muon local field H is proportional to the 
nearest moment size {H(z) oc S{z): local moment den- 
sity approximation), the distribution function p(H;Lq) 
of the local fields be obtained by the density of the states 
for the field H along the chain fragment Lq: 



P(H;L ) 



dz 



dH(z) 



(4) 



In the doped material, the impurity site should be lo- 
cated randomly. The experimentally observed field dis- 
tribution function p(H ) in this situation is the average of 
p(H; L ) for the Poisson distribution of the chain length 
Lq. The same procedure was taken for the model re- 
laxation employed in the analysis of the pSH spectra of 
doped CuGe0 3 [HI- 

The model muon relaxation function and the corre- 
sponding local field distribution function for the inhomo- 
geneous ordered moment size (Fig.QJ are shown in Fig.03 
The behavior of the distribution function p(H) is char- 
acterized by the ratio between the average chain-length 
A = (Lq) « 1/x and the moment-size recovery length £. 
In comparison to the model for doped CuGeC>3 (Fig. 4b 
of Ref. the local field distribution (Fig. [3Ji>) has the 

same shape but the field axis is reversed between H = 
and Hq. This is because in the doped CuGe03, im- 
purities create magnetic moments, whereas in Sr 2 CuC>3 
moments are depleted from the chain. The role of the 
zero-field H — and the maximum field H = Hq are re- 
versed between CuGeC>3 and Sr 2 Cu03 cases. As shown 
in Fig. [5^, when the ratio A/£ ^ 1, there is a coher- 
ent precession of muon spins because the moment size 
is relatively homogeneous and the local field distribution 
exhibits an isolated peak at H = Hq. As the ratio ap- 
proaches unity, the precession amplitude is diminished 
and at A/£ = 1, only a weak bump remains at the posi- 
tion of the first precession peak. Below A/£ = 0.5, there 
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FIG. 5: (a) The muon relaxation function and (b) local-field 
distribution for the inhomogeneously Neel ordered spin-chain 
(Fig. El . 



is no visible oscillation and the relaxation is mostly ex- 
ponential with time. Compared to the model for doped 
CuGe0 3 (Fig. 4b of Ref.jlj]) with the same level of back- 
ground relaxation, the amplitude of the precession signal 
component is much smaller in this model for depleted 
chains. The reason for the difference stems from the 
maximum moment size of the chain fragments: in the 
model of doped CuGe03, it was assumed that the ends 
of all chains have the same induced moment size sq. In 
this model for depleted chains, the maximum moment 
size appears at the center of the chain fragments, and 
depends on their length. 

If we represent the relaxation function shown in Fig.[S^, 
with a symbol P±(t; Hq, A/£), the muon relaxation func- 
tion for one muon site becomes: 

P M (t) = AjPxfatfo.A/O+Aiiexpft/Ti), (5) 

where A± {A\\) is the amplitude of the local field compo- 
nent which is perpendicular (parallel) to the initial muon 
spin polarization, and T\ is the spin-lattice relaxation 
time for the parallel component. As can be seen from the 
existence of two precession frequencies (Fig. |2K), there 
are two muon sites in Sr2Cu03 which, most likely, cor- 
respond to the muons attached to the in-chain and out- 
of-chain oxygen sites (see Appendix ^J. We introduced 
the local fields Hq and H B for the each site, and ob- 
tained the following phcnomcnological function for the 
muon spin relaxation: 



Ajf+ B ex P (VT 1 ), 



AfP^H^X/O 



(6) 



where A^_, A B and A^ +B arc the amplitudes of each 
components. Here, the two 1/T± signal components are 
combined in one relaxation amplitude (A^ +B ), since the 
background relaxation rate is too small to distinguish 
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Temperature (K) 

FIG. 6: (a) Muon spin relaxation in Sr2Cu03 analyzed with 
the relaxation function eq.© (solid line). The dot-dashed 
line is the contribution from the parallel- field term (A^ +B ) 
introduced in the analysis. Dashed lines are eq.© with the 
same parameters as in the x = case, except the A/£ ratio, 
(b) Temperature dependence of the recovery length £ and the 
magnetic order parameter as derived from the analysis. The 
dashed and the dot-dashed lines are the theoretical tempera- 
ture dependence of £ calculated from the experimental order 
parameter (see section ITTT1 . 



the contributions from the two sites. For the nominally 
pure sample, the precession signals can be analyzed with 
eq.©. The solid line in Fig. 01 is the result of the 
fit and the dot-dashed line is the contribution from the 
1/Ti relaxation term. It is noted that there are no ex- 
trinsic parameters introduced to describe the damping 
of the oscillation amplitude nor the early front-end re- 
laxation. However, eq.© based on the inhomogeneous 
moment size distribution (Fig. [3J describes the over-all 
feature of the muon spin relaxation fairly well. The ra- 
tio between the average chain-length A and the recovery 
length £ yields A/£ « 7.2 from the analysis. The fast 
front-end relaxation at t < lfis is also described by the 
perpendicular terms (Aj_ and A±). In this analysis, the 
1/Ti relaxation rate is not as larg e as the one obtained 
in the previous analysis of Ref . [lj] . 

The magnitude of the Curie term in the susceptibility 
for the nominally pure S^CuOa (Table QJ suggests the 
average chain-length for this nominally pure sample is 
A w 1/(2 x 0.04%) w 1000 lattices. This and the ratio 
A/£ = 7.2 obtained from the fitting analysis suggests that 



the recovery length of the moment size is £ w 150 lattices, 
which is at least one order of magnitude longer than in 
the Neel state of doped CuGe0 3 in which f sa 10 [2l|. 
The absence of a nominal spin-gap in the 5 = 1/2 anti- 
ferromagnetic chain is most likely the cause for the corre- 
lation length being longer than in the doped spin-Peierls 
compound CuGeOa. The temperature dependence of the 
recovery length (£) and the magnetic order parameter 
((s z ) oc precession frequency) are plotted in Fig. 03. The 
temperature dependence of £ is weak and possiblly ex- 
hibiting a slight increase at higher temperatures. 

Based on the magnitude of the Curie terms, one can es- 
timate the ratio A/£ for the Pd-doped compounds, which 
are summarized in Table[l]in parenthesis. The model re- 
laxation functions eq.@) for A/£ = 1 and 0.4 are shown 
as the dashed lines in Fig. These two parameter val- 
ues corresponds to the Pd 0.5% and 3% doped samples, 
respectively. Since the /xSR spectra of the doped systems 
do not exhibit the spectral features which were present 
in the x = sample, it is not possible to experimentally 
determine the parameters in eq.©. Here we assumed 
the same amplitudes (A^_, A B and A^ +B ), local fields 

(Hq and Hq) and 1/Ti relaxation rate as determined in 
the x = sample, and varied the A/£ parameter. This 
may be a good approximation in low doping in which the 
orientation of the local fields are approximately the same 
as in the nominally pure case. The calculated relaxation 
functions eq.|@ agree well with the behavior of the ex- 
perimentally obtained /iSR spectra of S^Cui-^Pd^Oa 
as shown in Fig. [SJi. Eventhough the /iSR spectra in 
the Pd-doped samples do not show precession, this might 
originate in the static Neel order with spatial distribution 
of moment-sizes. The decoupling measurement (Fig. |2Jd) 
confirms that the relaxation is static, which is a support 
to the idea that it originates from Neel order. Because of 
the long recovery length £, the small amount of impurity 
ions disturbs the ordered-moment size in a large area, 
so that the spin precession of the local magnetic probe 
becomes invisible. 

In the fj,SR time spectrum of the x = 0.5% sample, 
there is a hump in the measured polarization at ~ 1/xsec 
(Fig. HJi) , at the position of the first precession peak for 
the nominally pure compound. From this result, it is 
suggested that the maximum moment size Sq does not 
change in the sample with 0.5% Pd doping. For the 
x = 1.0% and 3.0% samples, such hump disappears, indi- 
cating that as the doping level increases, the parameter 
enters to the A/£ < 1 regime. The estimate of the A/£ 
value (Table P) confirms this conclusion. 

In Fig. [7| the Neel temperature of Sr2(Cu,Pd)03 and 
the theoretical curve are plotted. We estimate the 
average chain length L of Sr2(Cu,Pd)03 from the Curie 
term of the susceptibility, assuming that one Pd ion cre- 
ates either S = or S = 1/2 with the equal probability. 
The horizontal and vertical axes of Fig. [7] are normalized 
by the inter-chain interaction J' which has not been ex- 
perimentally obtained for Sr2Cu03. However the ratio 
Tn/J = 4 x 10~ 4 has been estimated an d the tue ~ 
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FIG. 7: Neel temperature vs. average chain length. The 
theoretical line is taken from Eggert et al. [l^l . 



ory proposes Tjy / J' ~ 2.6 in the long chain-length limit 
(Fig. 0) . From these, one can obtain the ratio between 
the inter- and intra-chain couplings J'/ J » 1.6 x 10~ 4 . 
This value was employed to scale the average chain length 
L in the horizontal axis of Fig. [7] The agreement between 
theory and measurement regarding the Neel temperature 
is reasonablly good, suggesting that the theory correctly 
estimates the energy scale of the magnetic order. 

Within the same theoretical framework, the spatial dis- 
tribution of the ordered moment size has been calculated 
recently j^. The assumption made is that the S = 1/2 
antiferromagnetic chain acquires the staggered moment 
size (s z ) in a consistent way from the effective staggered 
magnetic field B originating from the neighboring or- 
dered chains. The intra-chain interaction J propagates 
the effect of the chain-end, reducing the ordered moment 
size at its vicinity. The inter-chain interaction J' recov- 
ers the ordered moment size by inducing the staggered 
magnetic field from the neighbouring chains. The com- 
petition between these two interactions on the chain frag- 
ments causes the spatial distribution of ordered moment 
sizes as shown in Fig. QJ A scaling argument proposes 
that the recovery length £ is given by £ oc i?~ 2 / 3 |22ll23|. 
where B is the effective staggered field. In the mean- 
field approximation, the staggered field can be written 
as B = zJ'(s z ), where J' is the inter-chain interaction, z 
is the number of neighboring chains and (s z ) is the mag- 
netic order parameter. Combining these two relations, 
the recovery length in the mean-field approximation is 
given by £mf k (s z )~ 2 ^ 3 - There is also a scaling relation 
between the magnetic order parameter and the staggered 
field: (s x ) oc B 1 / 3 22]. With this and the scaling rela- 
tions relation between £ and B, one can obtain a scaling 
prediction for the recovery length: £sc cx (s z )~ 2 - In 
either case, the temperature dependence of the recovery 
length is set by that of the order parameter (s z ) but with 



a different exponent. 

The temperature dependence of the recovery length in 
the mean-field approximation (£mf) and the scaling re- 
sult (£sc) are shown in Fig. Et>, based on the experimen- 
tally obtained magnetic order parameter. Since we do not 
know the over-all scale factor of the recovery length, we 
employed the experimental value obtained at low temper- 
ature. The experimental result is more consistent with 
the full scaling result £sc, rather than with the value ex- 
pected in the mean-field framework, except at the vicinity 
of the Neel temperature. The critical exponent of the or- 
der parameter as a function of temperature (/3 rj 0.2) [2^ 
also indicates that the system is in the scaling regime 
rather than in the three-dimensional mean-field regime 
in which (3 w 0.5. Our measurements suggest that the 
finite recovery length of the ordered moment sizes is de- 
termined by the effective staggered magnetic field B orig- 
inating from the neighboring chains. The over-all scale 
factor £/a « 150 in the low temperature limit is yet to be 
calculated theoretically for the parameters of S^CuOs. 

In the theoretical calculation of chain fragments with 
the lengths longer than the recovery length (Lo/£ ft l)j 
the self-consistent moment size exhibits a two peak struc- 
ture in the distribution function (Fig. 1 of Ref. 22]). One 
peak is located almost at zero-moment size and the other 
appears at a slightly larger size than that for the unper- 
turbed chain. The former and the latter peaks originate 
from the chain fra gme nts with even and odd number of 
spins, respectively [22]. We have found that the Fourier 
transformation of the moment size distribution function 
exhibits clear precession, as a result of the isolated peak 
for the odd-numbered fragments. This theoretical result 
is in contrast to the experimental observation (Fig. [2K) , 
in which the precession is dampled upon doping. There 
are two possible explanations: (1) The moment size of 
the odd length chains may be reduced due to quantum 
fluctuations, as has already been suggested in Ref. [2^] . 
(2) The sum of dipolar fields at the muon site contains 
contributions from neighboring chains, smearing out the 
even-odd effect. The proposed muon sites and the dipo- 
lar field calculations ( Appendix |A"|) demonstrate that the 
contribution from the nearest neighbor chain is actually 
dominant, at least for the higher frequency site. This 
suggests that the quantum fluctuation scenario (1) may 
be a more favoured explanation for the absence of muon 
spin precession in Sr2(Cu,Pd)03. 



IV. CONCLUSIONS 

We have investigated how Neel order is destructed 
with non-magnetic impurity doping in the quasi one- 
dimensional S — 1/2 antifcrromagnet. The model ma- 
terial employed is Sr 2 (Cui_ x Pd K )03. The susceptibility 
at T > Tm exhibits Curie- Weiss behavior. The mag- 
nitude of the Curie term is consistent with a model in 
which the creation of half an induced moment (size = 
1/is) occurs due to substitution with one Pd ion. This 



suggests that the induced paramagnetic moment origi- 
nates from the total spin of the chain fragments which 
is either S — or S — 1/2 depending on the fragment 
consists of even or odd number of spin sites. The muon 
spin relaxation of the nominally pure Sr2Cu03 was re- 
analyzed with a model for the spatially inhomogeneous 
Neel ordered state. The length scale £ was introduced 
to describe the recovery of ordered moment size away 
from the spin defect. It was suggested that the length 
scale £ originates from the effective staggered magnetic 
field of the neighboring chains. The length scale £ w 150 
lattices for S^CuOa is more than 10 times longer than 
in doped CuGeC>3, which reflects the nominally gapless 
characteristic of the system. 

In the Pd doped samples, the absence of muon spin pre- 
cession for less than 1% doping level is consistent with 
a large recovery length £. The Neel temperature, which 
was defined as the temperature where muon spin relax- 
ation rate starts to increase, exhibits a good agreement 
with recent theoretical calculation. 
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APPENDIX A: MUON SITES 

The crystal structure of Sr2Cu03 has two non- 
equivalent oxygen sites [l^. Site 0(1) is out of the chain 
oxygen site of the corner-shared CUO3 plaquett, and site 
0(2) is the shared oxygen of the two neighboring pla- 
quetts forming the chain structure. In cuprates, it has 
been proposed that muons form an 0-/i + bond with oxy- 
gen ions in an analogy with the hydrogen bonding |24j . 
In this appendix, we calculate the electro-static potential 
for a muon which was assumed to form an O-fj, bond with 
a bond length of lA and determine the orientation of the 
bond with respect to the crystallographic axis. We also 
calculate the dipolar fields and compare with the exper- 
imentally observed local fields. 

In Fig.[Hk L andlHb, contour plots of the electro-static po- 
tential are drawn for the muons forming the O-fi bonds 
with the 0(1) site and 0(2) site, respectively. The 
electro-static potential is calculated using the method of 
Ewald sum |2f| assuming the formal point charges lo- 
cated at the ionic positions. The minimum of the po- 
tential exists at (8,ip)^ = (60°, 90°) and (0, <p)(2) = 
(20°, 0°) for the 0(1)-// and 0(2)-/z bonds, respectively, 
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O-ji bond orientation if) (deg) O-ji bond orientation if) (deg) 



FIG. 8: (a) and (b): Electro-static potential for a muon lA 
away from the oxygen site O(l) and 0(2), respectively. The 
unit of the potential is in eV. The crosses are the potential 
minima which are the most probable orientation for the O- 
M bond, (c) and (d): Dipolar fields calculated for the same 
O-/X bonding. The unit of the field is in Gauss for the spin 
structure determined by the neutron diffraction measurement 
0- 



which are shown by the cross symbol in the figures. The 
polar coordinate is defined as the = direction be- 
ing parallel to the longest crystallographic axis (c-axis), 
and the (p = 0° direction in the basal plane being par- 
allel to the CuO chains (a-axis). In Fig. the muon 
sites in the real space crystal structure are shown by the 
star-symbols. The potential minimum value for the bond 
length of lA is deeper for the 0(l)-/i site (-10.4eV), than 
for the 0(2)-/i site (-9.2eV). This suggests that 0(2)-^ 
site may have a shorter bond-length to gain the electro- 
static potential of the oxygen ion. 

The ordered moment size and its orientatin of Sr2Cu03 
has been obtained by neutron diffraction measurements 
as 0.06/is pointing along the CuO chain direction, re- 
spectively [14|. The corresponding dipolar fields for the 
spin structure obtained are calculated for the two kinds 
of 0-/J, bonds. The results are shown in Fig. [Ht and|HJi. 
At the orientation of the minimum potential (cross sym- 
bols in the figures), the dipolar fields are 24G and 150G 
for the 0(l)-/z and 0(2)-/Lt site, respectively. The ex- 
perimentally observed local fields in S^CuOs are 23. 3G 
and 97. 7G [14|. The former is close to the calculation of 
0(l)-^i site, which is most likely the muon position re- 
sponsible to the 23. 3G signal. The 97. 7G signal does not 
agree well with the calculation for the 0(2)-/j, site. The 
disrepancy might originate from the bond length which 
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FIG. 9: (a) The crystal structure of Sr 2 Cu0 3 . 0(l)-fi and 
0(2)-/i sites as determined by the electro-static potential cal- 
culation are shown by the star symbols which are attached to 
the each oxygen sites, (b) The dipolar fields from the each 
chains to the muon sites. The boxes are distributed at the 
position of the chains seen along the a-axis. The upper and 
lower panels of the boxes are the dipolar fields for 0(l)-/t and 
0(2)-fi sites, respectively. The unit of the fields is in Gauss. 



could be shorter for 0(2)-/i, site. We have calculated 
the electro-static potential and the dipolar fields for a 
shortened 0(2)-/i bond length of 0.9A and found that 
the dipolar field at the potential minimum is reduced to 
120G for the same spin structure. Since the pseudo min- 
imum position for the 0(l)-/z site (60°,0°)( 1 ) exhibits a 
much smaller dipolar field (40G), the 97. 7G signal most 
likely originates from the 0(2)-/i site with a shortened 
bond length. 



For the muon sites proposed above, we calculate the 
contribution of the dipolar fields from neighboring chains. 
In Fig. I^Ja, we show the magnitude of the dipolar fields 
from the chains at the position of the boxes which rep- 
resent the CuO chains projected to the b x c-plane. The 
upper and lower panels of the boxes indicates the dipolar 
fields of the 0(l)-/x and 0(2)-/x sites, respectively. For 
the 0(l)-/x site, the contribution from the nearest neigh- 
boring chain is as large as 1/3 of the main contribution. 
This contribution from the neighboring chains broadens 
the local field distribution from that expected by the spa- 
tial distribution of the ordered moment sizes. However, 
the precession signal with higher frequency of the nom- 
inally pure Sr2Cu03 originates from the 0(2)-// site, at 
which the contribution of the dipolar field from the neigh- 
boring chains is at most 1/10 of the main contribution. 
This calculation suggests that the precession signal with 
higher frequency would survive, if the ordered moment 
size distribution has a distictive peak at a finite frequency 
as proposed in the theoretical calculation j^. The ab- 
sense of the precession signal in the Pd-doped sample 
suggests that the isolated peak of the ordered moment 
size, which originates from the chain fragments with odd 
number of spins, is actually diminished in Sr2Cu03 due 
to the effects which are not included in the theoretical 
calculation. 
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